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Evolving faults are starting in one phase of circuit and spreading to other 
phases after some time. There has not been a suitable method for locating 
evolving faults in double circuit transmission line until now. In this paper, 
a novel method for locating different types of evolving faults occurring in 
double circuit transmission line is proposed by considering adaptive 


neuro-fuzzy inference system. The fundamental current and voltage 


magnitudes are specified as inputs to the proposed method. The simulation 
results using MATLAB verify the effectiveness and correctness of the 
protection method. Simulation results show the robustness of the method 
against different fault locations, resistances, time intervals, and all evolving 
fault types. Moreover, the proposed method yields satisfactory performance 
against percentage errors and fault location line parameters. The proposed 
method is easy to implement and cost-effective for new and existing double 
circuit transmission line installation. 
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1. INTRODUCTION 

The growth of a country depends on availability of continuous electrical power supply. 
The power system consists of generating, transporting and utilizing electrical power. The only means of 
electrical power transmission is through transmission lines. Most of the transmission lines are overhead 
transmission lines which are open to the surrounding. Transmission lines are subject to different faults 
due to snow, storms, insulators breakdown, lightning, large birds and other external objects. 
Besides, when shunt fault takes place on transmission line in rural places, bushfire hazards may occur. 
Thus, detection of fault location within transmission line is crucial to minimize revenue loss, 
reduce unplanned outage time, improve network reliability and restore service rapidly. Many transmission 
line fault locating schemes based on soft computing techniques and double circuit transmission analysis have 
been introduced in the following literature work. 

Most of the research work has been carried out for protection of double circuit three phase 
transmission lines (DCTPTL) [1-5]. The location of faults by considering fundamental component 
of currents and voltages in transmission line has been widely studied in Ref [6]. Further, different fault 
location algorithms have been presented in DCTPTL based on fuzzy inference system (FIS) [7-8], 
support vector machines [9], artificial neural networks (ANN) [10] and k-nearest neighbor algorithm [11]. 
Presence of a variety of individual computing methods [7-11] available for locating the DCTPTL is still 
having certain problems which limit their applicability. Hybridization of technologies can solve the 
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limitations of individual computing applications. The adaptive neuro-fuzzy inference system (ANFIS) can be 
characterized as a hybrid’ method, as it represents the integration of ANN and fuzzy logic technique. 
Recently, ANFIS has been found highly successful in applications such as evaluation of drought indices [12] 
and prediction of solar thermal energy system [13]. 
In the last few decades, ANFIS networks have been extensively used by researchers for location of 
transmission line faults [14-20]. Although, all the ANFIS location schemes [14-20] are designed to line to 
line faults (normal shunt faults) occurring at single time but are not applicable for line to line faults 
(evolving faults) occurring with time gap. Evolving faults include two sequential faults, resulting in the 
change of fault phase and time in DCTPTL. Some notable works have been reported in the literature on the 
location of evolving faults in transmission lines [21-27]. The complexity arising due to the time gap of two 
faults which adversely affect the performance of transmission line relaying. However, there is still a need to 
design fault location method to handle evolving faults. Hence, it can be concluded that, hitherto none of the 
earlier reported papers can locate evolving faults using ANFIS and DCTPTL current and voltage magnitudes. 
In this regard, a novel method to locate evolving faults in the DCTPTL using ANFIS is presented 
in this paper. The novelty and highlights of the proposed work can be summarized as 
1. Development of ANFIS based fault locator using single end fundamental current and voltage 
magnitudes and thus avoids the use of communication channel for sending and receiving terminals 

2. Development of ANFIS without using classification method 

3. Detects the evolving faults correctly 

4. Further, another notable contribution of the proposed method is that it improves the accuracy of fault 
location and 

5. Effectively distinguish the two successive faults different time positions on transmission lines. 

The rest of this paper is structured as follows: Section 2 describes the DCTPTL simulation and 
training pattern generation. Section 3 presents the fault location algorithm for evolving faults using ANFIS. 
Section 4 reports the simulation results of fault location method. Section 5 provides conclusions. 


2. TRANSMISSION LINE SIMULATION AND TRAINING PATTERN GENERATION 

Figure 1 depicted the simplified diagram of the chosen 100 km, 400 kV, 50 Hz, DCTPTL. 
The MATLAB software is employed for DCTPTL and evolving faults simulation. DCTPTL is constructed 
using distributed parameter line block. Different evolving faults considered in this paper are | phase to 
ground fault changed to 2 phase to ground, 2 phase to ground faults changed to 3 phase to ground, 
2 phase fault changed to 3 phase to ground fault and lastly 1 phase to ground fault changed to 2 phase to 
ground which again changed to 3 phase to ground. 


Bus-1 Bus-2 


—— 100 km ——> 


Source-S : . Source-R 
400 kv Circuit-2 400 kV 
I Load-2 


Figure 1 The simplified diagram of DCTPTL 


The three phase current in circuit-1 JA1, IB1, IC1, and voltage VA, VB, VC waveforms of corresponding 
phases during an evolving fault in phase ‘A’ at 72 ms time and phase ‘B’ at 84 ms time at 80 km from the 
bus-1 are exemplified in Figure 2 and Figure 3. 
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Figure 2. Current signals obtained during evolving faults 
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Figure 3. Voltage signals obtained during evolving faults 


3. PROPOSED METHOD 

The fuzzy and ANN have their own benefits and drawbacks. ANFIS has capability to capture the 
merits of both in a single construction. It is a type of adaptive network that is equivalent to FIS functionally. 
In this, to improve the accuracy of FIS model, neural network is used. The benefits of ANFIS are excellent 
explanation facilities through IF-Then rules, fast, easy to implement, accurate learning, easy to integrate both 
linguistic and numeric knowledge and strong generalization abilities. ANFIS can be adopted in a variety of 
applications of designing, signal processing, control and decision making. It can also be employed for both 
nonlinear and linear of input-output parameters. 

For training of the ANFIS fault locator, various fault cases are simulated with change in evolving 
fault type, time interval, fault location and fault resistance. Three-phase voltages and currents in circuit-1 and 
circuit-2 are collected from the source point for all fault cases. The fault parameters used for training of the 
ANHIS are listed in Table 1. It can be observed from the Table | that the number of fault cases for training is 
1728. After simulating various fault cases with changing fault parameters, the voltage and current magnitudes 
measured at the source point are processed. 

The ANFIS has constructed in MATLAB software with the DCTPTL simulation data. 
A single ANFIS has been developed to locate all types of evolving faults in both the circuits under different 
DCTPTL operating conditions. Nine inputs and single output are taken for developing ANFIS. 
The fundamental component of three phase current and voltage amplitudes are fed to ANFIS for locating the 
evolving faults. The fuzzy system inputs symbolize IAIf, IB1f, ICIf, IA2f, IB2f, IC2f, VAf, VBf and VCf. 
Its output symbolizes location of fault Lf. All input current values are classified as Il to IS. 
Similarly, each voltage inputs are classified as V1 to V5 and output is classified as Lfl to Lf5. The inputs and 
outputs are classified based on Gaussian membership functions. An array of five IF-Then rules is made for 
fuzzy system as explained below. 


IF-THEN Rules 

Rule-1:IF (IA1f is [5) or (IBIf is I5) or (ICIf is I5) or (IA2f is 15) or (IB2f is I5) or (IC2f is 15) or (VAF is 
V1) or (VBf is V1) or (VCf is V1) THEN (Lf is Lf1) 

Rule-2:IF (IA1f is 14) or (IBIf is I4) or (ICIf is 14) or (IA2f is 14) or (IB2f is 14) or (IC2f is 14) or (VAF is 
V2) or (VBf is V2) or (VCf is V2) THEN (Lf is Lf2) 
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Rule-3:IF (IA1f is 13) or (IBIf is I3) or (ICIf is I3) or (IA2f is 13) or (IB2f is 13) or (IC2f is 13) or (VAF is 
V3) or (VBf is V3) or (VCf is V3) THEN (Lf is Lf3) 
Rule-4:IF (IA1f is I2) or (IBIf is I2) or (ICIf is I2) or (IA2f is I2) or (IB2f is I2) or (IC2f is 12) or (VAF is 
V4) or (VBf is V4) or (VCf is V4) THEN (Lf is Lf4) 
Rule-5:IF (IA1f is I1) or (IB1f is 11) or (ICIf is I1) or (1A2f is 11) or (IB2f is I1) or (IC2f is 11) or (VAF is 
V5) or (VBf is V5) or (VCf is V5) THEN (Lf is Lf5) 


The fuzzy inference is computed by Takagi Sugeno model. The FIS and training data are given to the ANFIS. 
The ANN is training by back-propagation method. The most excellent performance is achieved by using 
back-propagation algorithm with 80 epochs, 5 Gaussian membership functions in each input and output with 
5 IF then rules has been used. The structure of the developed ANFIS model is shown in Figure 4. In this way, 
ANFIS is developed, which provides only one output for location of evolving faults. Now ANFIS is ready to 
collect the current and voltage and can find the location of evolving faults. When the current and voltage 
magnitudes are given to the developed method, one of the aforementioned rules is fired and the 
corresponding location of evolving faults is obtained. 


Input Inputmf Rule Outputmf Output 
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Figure 4. The structure of the developed ANFIS 


Table 1. Training parameter changing 
Training data 


Fault parameters 


A1-G to AIB1-G, 
B1-G to BIC1-G, 
C1-G toC1A1-G, 
A2-G to A2B2-G, 
B2-G to B2C2-G, 
C2-G to C2A2-G. 
AIB1 to AIBAIC1-G, 
BIC1 to AIBAICI1-G, 
C1A1 to AIBAIC1-G, 


1 phase to ground fault changed to 2 phase to ground fault 


2 phase fault changed to 3 phase to ground fault 


Fault type 


A2B2 to A2BA2C2-G, 
B2C2 to A2BA2C2-G, 
C2A2 to A2BA2C2-G. 


(Evolving faults) 
[24] 


2 phase to ground faults changed to 3 phase to ground fault 


A1B1-G to AIBAIC1-G, 
BIC1-G to AIBAIC1-G, 
C1A1-G to AIBAIC1-G, 
A2B2-G to A2BA2C2-G, 


B2C2-G to A2BA2C2-G, 
C2A2-G to A2BA2C2-G. 
A1-G to A1B1-G to AIBIC1-G, 
B1-G to BIC1-G to AIBI1C1-G, 
C1-G to C1A1-G to AIBICI1-G, 
A2-G to A2B2-G to A2B2C2-G, 
B2-G to B2C2-G to A2B2C2-G, 
C2-G to C2A2-G to A2B2C2-G. 
Fault location (km) [12] 1,5, 15, 25, 35, 45, 55, 65, 75, 85, 95, 99. 

Fault resistance (Q) [2] 0, 100. 
Time interval (ms) [3] 10, 20, 30. 
Total number cases 24x12x2x3= 1728 


1 phase to ground fault changed to 2 phase to ground which 
again changed to 3 phase to ground fault 
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4. RESULTS AND DISCUSSION 

The performance of the proposed ANFIS based fault location method has been examined and the 
results have been analyzed. In order to examine the response of proposed approach, it is necessary to get the 
relevant data which is never used during training process. Different parameters have been selected for 
checking the proposed ANFIS. Total 14000 test samples of evolving faults, with change in fault parameters 
like faults locations, fault resistances, fault types and time intervals have been collected. After the training 
process, the achieved optimal network configurations are checked with the MATLAB software tool samples. 
The error (%) in estimated fault location is calculated using (1). Some results of the evolving faults and 
proposed method are explained below. 


Esitimated fault location (km)—Real fault location (km) 


Error (% ) = | | x 100 (1) 


Line length 


4.1. Performance changing evolving fault types 

In addition to ten types of shunts faults, DCTPTL is also prone to evolving faults due to 
animal/human and tree branch falling on live conductor normally begins off as single phase to ground faults. 
Most of the methods have been developed and unable to located especially evolving faults. 
Thus, the developed method is verified for different evolving faults in both the circuits of transmission lines. 
Performance changing evolving fault types, keeping invariable time interval (20 ms), constant 
resistance (80 Q), and fixed location of fault (65 km) are listed in Table 2. The error is up to 0.281% and 
detection time is up to 0.06 ms for all types of evolving fault. As seen in Table 2, evolving faults have no 
influence on the performance of the proposed method even for different evolving faults occurring in the line. 


Table 2 Performance changing evolving fault types 


Initial Evolving Time Fault Resistance — Fault Location Estimated Error Detection Time 
Fault Fault Interval (ms) (Q) (km) Location (km) (%) (ms) 
Al-G A1B1-G 20 80 65 64.940 0.060 0.02 
B1-G BIC1-G 20 80 65 65.252 0.252 0.03 
C1-G C1A1-G 20 80 65 65.082 0.082 0.01 
A1B1 AIBAICI1-G 20 80 65 65.201 0.201 0.06 
BIC1 AIBAIC1-G 20 80 65 64.877 0.123 0.01 
C1AI AIBAIC1-G 20 80 65 65.011 0.011 0.04 
AIBI-G  AIBAICI1-G 20 80 65 65.143 0.143 0.06 
BIC1-G = AIBAIC1-G 20 80 65 64.841 0.159 0.02 
C1AI-G = AIBAICI-G 20 80 65 65.160 0.160 0.04 
Al1-G AIBICI-G 20 80 65 65.105 0.105 0.06 
B1-G AIBICI-G 20 80 65 64.841 0.159 0.04 
C1-G AIBICI-G 20 80 65 65.063 0.063 0.04 
A2-G A2B2-G 20 80 65 65.212 0.212 0.02 
B2-G B2C2-G 20 80 65 64.903 0.097 0.01 
C2-G C2A2-G 20 80 65 64.021 0.021 0.03 
A2B2 A2BA2C2-G 20 80 65 65.222 0.222 0.01 
B2C2 A2BA2C2-G 20 80 65 65.029 0.029 0.02 
C2A2 A2BA2C2-G 20 80 65 65.098 0.098 0.06 
A2B2-G  A2BA2C2-G 20 80 65 65.281 0.281 0.04 
B2C2-G  A2BA2C2-G 20 80 65 65.025 0.085 0.02 
C2A2-G  A2BA2C2-G 20 80 65 65.192 0.192 0.06 
A2-G A2B2C2-G 20 80 65 65.144 0.144 0.02 
B2-G A2B2C2-G 20 80 65 64.942 0.058 0.03 
C2-G A2B2C2-G 20 80 65 65.169 0.169 0.02 


4.2. Performance changing time intervals 

The faults can take place with evolving fault type at any time; therefore, to study its effect on the 
developed method, a variety of time intervals for evolving faults have been applied. In order to explain the 
impact of time interval, the fault in the fixed location value (43km) with constant resistance value (50Q) for 
the evolving fault type (A2-G to A2B2-G) and various time intervals has been considered. Some of results of 
fault location changing time intervals during evolving faults are illustrated in Table 3. The error is up to 
0.290% and detection time is up to 0.07 ms for all cases. The maximum error of the fault location in 
DCTPTL line remains acceptable regardless of the variation in the time intervals. 
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Table 3 Performance changing time intervals 


Initial Evolving Time Fault Fault Estimated Error (%) Detection Time 
Fault Fault Interval (ms) Resistance (Q) Location (km) Location (km) (ms) 
A2-G A2B2-G 21 50 43 42.845 0.155 0.04 
A2-G A2B2-G 26 50 43 43.122 0.122 0.03 
A2-G A2B2-G 19 50 43 43.253 0.253 0.01 
A2-G A2B2-G 10 50 43 43.219 0.219 0.02 
A2-G A2B2-G 17 50 43 43.290 0.290 0.02 
A2-G A2B2-G 15 50 43 42.948 0.052 0.01 
A2-G A2B2-G 28 50 43 42.897 0.103 0.07 
A2-G A2B2-G 18 50 43 43.186 0.186 0.05 
A2-G A2B2-G 16 50 43 43.096 0.096 0.05 
A2-G A2B2-G 25 50 43 42.759 0.241 0.01 
A2-G A2B2-G 17 50 43 43.010 0.010 0.03 
A2-G A2B2-G 11 50 43 43.120 0.120 0.03 


4.3. Performance changing fault resistance 

The fault resistance has a significant effect on the study of evolving fault analysis. Various fault 
resistances have been used in the simulation study and the results with typical fault resistances are reported 
here. To check its performance the evolving fault (AIB1-G to AIBA1C1-G), kept constant location (18km) 
and time interval (15ms), and resistances are varied from 10Q to 120Q. The test results of the developed 
method for evolving faults in DCTPTL under changing resistances are illustrated in Table 4. The error is up 
to 0.230% and detection time is up to 0.06 ms for all types of fault. Results give a satisfactory error of the 
proposed method for different fault resistances. What’s more, the method stays unaffected by fault resistance. 
Even if the fault occurs on the vicinity of resistances, the method can locate evolving faults. 


Table 4. Performance changing fault resistance 


Initial Evolving Fault Time Fault Resistance Fault Estimated Error Detection Time 

Fault Interval (ms) (Q) Location (km) Location (km) (%) (ms) 
A1B1-G AIBAIC1-G 15 5 18 18.078 0.078 0.03 
AI1B1-G AIBAIC1-G 15 10 18 18.056 0.056 0.04 
AI1B1-G AIBAIC1-G 15 20 18 17.770 0.230 0.05 
AI1B1-G AIBAIC1-G 15 20 18 18.198 0.198 0.03 
AI1B1-G AIBAIC1-G 15 40 18 18.206 0.206 0.02 
AI1B1-G AIBAIC1-G 15 50 18 18.135 0.135 0.02 
AI1B1-G AIBAICI1-G 15 60 18 17.821 0.179 0.05 
AI1B1-G AIBAIC1-G 15 70 18 17.838 0.162 0.03 
AI1B1-G AIBAIC1-G 15 80 18 18.047 0.047 0.04 
A1B1-G AIBAIC1-G 15 90 18 18.228 0.228 0.06 
AI1B1-G AIBAIC1-G 15 95 18 17.946 0.054 0.01 
AI1B1-G AIBAICI-G 15 99 18 18.101 0.101 0.02 


4.4. Performance changing fault location 

In this part, the simulation results with different fault locations are presented. The developed method 
is tested by keeping fixed time interval (10ms), constant resistances (20Q), different locations of evolving 
fault (B1-G to B1C1-G) from the relaying point ranging from Lf = 1 km to 99 km of the line. Some of the 
results of evolving faults in DCTPTL under changing locations are given in Table 5. The maximum absolute 
error percentage was 0.292% under evolving faults and maximum detection time taken by the proposed 
method is 0.06 ms. As can be seen in Tables 5 the developed method is able to achieve a reasonably high 
accuracy in locating fault. Also, the performance of the developed method is not affected by variations in the 
fault location. 
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Table 5 Performance changing fault location 

Initial Evolving Time Fault Resistance Fault Estimated Error Detection 
Fault Fault Interval (ms) (Q) Location (km) Location (km) (%) Time (ms) 
B1-G BICI-G 10 20 9 8.926 0.074 0.05 
B1-G BICI-G 10 20 15 15.085 0.085 0.03 
B1-G BICI-G 10 20 28 27.068 0.068 0.01 
B1-G BICI-G 10 20 36 35.891 0.009 0.01 
B1-G BICI-G 10 20 44 44.104 0.104 0.02 
B1-G BICI-G 10 20 51 51.211 0.211 0.03 
B1-G BICI-G 10 20 60 60.282 0.282 0.02 
B1-G BICI-G 10 20 67 66.135 0.135 0.01 
B1-G BICI-G 10 20 73 72.823 0.177 0.06 
B1-G BICI-G 10 20 82 82.292 0.292 0.03 
B1-G BICI-G 10 20 89 89.124 0.124 0.05 
B1-G BICI-G 10 20 99 98.761 0.239 0.04 

5. CONCLUSION 


This paper has outlined a new ANFIS method for locating evolving faults as well as shunt faults in 


the DCTPTL. System designing and simulations have performed in MATLAB software package for 
capturing samples of faulty signals by varying fault parameters. Low pass filter, sampler, and DFT have 
utilized to derive the features of raw current and voltage signals, which are further given as inputs for the 
fault locator. The proposed technique has overcome problems such as information from two ends of the 
transmission line and knowing the fault type. The simulation results showed that the method can be applied 
within the entire line and the accuracy of the fault location is high. Furthermore, the proposed method 
correctly located the evolving faults and tested with all fault events. Thus it can be applied for protection of 
real power systems as well. 
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